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Abstract 

We obtain an improved upper bound on the lightest neutral CP-even Higgs boson 
mass from a low energy renormalisation-group analysis of the Higgs sector of the next- 
to-minimal supersymmetric standard model. We find rrih < 145 GeV for m t = 90 
GeV, decreasing to m h < 123 GeV for m t = 180 GeV. We also discuss the light Higgs 
spectrum in the region of parameter space close to the upper bound. 



The question of the origin of electroweak mass is one of the most urgent ques- 
tions of present day particle physics. The discovery of a particle which resembles 
the Higgs boson of the minimal standard model, and the measurement of its mass, 
will provide clues as to the nature of new physics beyond the standard model. In 
this letter we shall be concerned with the question of how heavy the lightest neutral 
CP-even supersymmetric Higgs boson, h°, can be within the framework of super- 
symmetric grand unified theories (SUSY GUTs) Q. In SUSY GUTs all the Yukawa 
couplings are constrained to remain perturbative in the region Msusy ~ 1 TeV to 
Mqjjt ~ 10 16 GeV. This constraint provides a maximum value at low energies for 
those Yukawa couplings which are not asymptotically-free, and is obtained from the 
renormalisation group (RG) equations together with the boundary conditions that 
the couplings become non-perturbative at Mqut ~~ the so-called "triviality limit" . In 
the minimal supersymmetric standard model (MSSM) M, the triviality limits provide 
a useful bound on the top quark mass m t . The upper bound on the h° mass, m h , 
in the MSSM, including radiative corrections, has recently been the subject of much 
discussion^], |], ||, |7], M. However the MSSM is not the most general low energy 
manifestation of SUSY GUTs. It is possible that SUSY GUTs give rise to a low energy 
theory which contains an additional gauge singlet field, the so called next-to-minimal 



supersymmetric standard model (NMSSM) [|H], [Ll], 0. Here we shall concentrate 
on the question of the upper-bound on rrih in the NMSSM which is obtained from 
triviality limits of Yukawa couplings. 

In the NMSSM there are two Higgs doublets Hi, H 2 and a complex Higgs singlet 
N, leading to three neutral CP-even Higgs bosons, two neutral CP-odd states and 
two charged bosons in the physical spectrum. The superpotential has the form 

W = h t Qt c H 2 + XNHiH 2 - ^N 3 + . . . , (1) 

where the superfield Q contains the left-handed third family quark doublet 
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the superfield t c contains the charge conjugate of the right-handed top quark field, 
H\ 2 contain the Higgs doublets and N contains the Higgs gauge singlet. The el- 
lipsis represents terms whose relatively small couplings will not play a role in our 
analysis. The NMSSM removes the fj,HiH 2 term in the superpotential of the MSSM 
and replaces its effect by the vacuum expectation value (VEV) of the Higgs singlet, 



< iV >= x. The other VEVs are < Hi 2 >= 1*1,2, where v = \fv\ 2 + v 2 2 = 174 GeV. 
The trilinear term |iV 3 is necessary in order to avoid an approximate global U(l) 
symmetry (broken by instanton effects) which, when spontaneously broken by the 
VEVs, would lead to an axion. The Higgs sector of the NMSSM reduces to that of 
the MSSM in the limit r = x/v — > 00 with Ax and kx fixed, where one of the neutral 



CP-even states and one of the CP-odd states decouple []TT . 

An upper bound on the lightest neutral CP-even scalar h° in the NMSSM may 
be obtained from the real symmetric 3x3 neutral scalar mass squared matrix, by 
using the fact that the smallest eigenvalue of such a matrix must be smaller than the 



smallest eigenvalue of its upper 2x2 block. The resulting bound at tree-level is |12 



mh < M z 2 + (A V - M z 2 ) sin 2 2(3. (2) 

where tan/3 = and A is regarded as a running parameter evaluated at M S usy- The 
upper bound on is determined by the maximum value of X(Msusy), henceforth 
denoted X max . The value of X max is obtained by solving the SUSY RG equations for 
the Yukawa couplings h t , A and k in the region Msusy = 1 TeV to Mqut = 10 16 GeV 
fl~3| , PI] . If the Yukawa couplings h t , A and k are all initially large at Mqut then they 
approach low energy fixed point ratios However, if the boundary condition at 

Mqut is A ^> k, then larger values of X(Msusy) can be achieved. We have repeated 
the calculation of ref.[l4| and found that for h t (M SUSY ) = 0.5 — 1.0, X max = 0.87—0.70 



and for h t (MsusY) 1.06, Q X max — > (with k = always). Radiative corrections to 

1 h t (MsuSY) < 1-06 is the triviality bound, which, together with m t — ht(mt)v sin /3 ', where 
ht{m t ) < 1.12, implies the bound m t < 195 GeV. 



the tree-level bound in Eq.(2) have been considered in refs. ||T5|, |I6, |I7|]. In ref.|l6[ these 



were estimated from a low energy RG analysis of the Higgs sector of the model between 
Msusy and a lower scale /i, assuming that only one Higgs boson has a mass below 
Msusy- Here we shall consider the more general case in which both Higgs doublets 
and the Higgs singlet may be lighter than Msusy, making the usual approximation 
of hard decoupling below M S usy of the superpartners. We shall then apply this 
technique to obtain a new upper bound on rrih, and supplement the discussion with 
a full numerical analysis of the light Higgs spectrum of the NMSSM in the region of 
parameter space close to this upper bound. 

At energy scales above Msusy the Higgs potential in the NMSSM with general 
soft SUSY breaking terms is given from Eq.(l) by 
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H ± \ 2 + \H 2 \ 2 )\N\ 2 + \H X H 2 \ 2 + £; 2 |iV| 4 - Xk(HiH 2 N* 2 + H.c.) 



+ l -{ 9l 2 + g 2 2 ){\H 2 \ 2 - IH,] 2 ) 2 + 

+ ttih 2 ^ 2 + m H 2 \H 2 \ 2 + m N 2 \N\ 2 

- \A\(HiH 2 N + H.c.) - hA k (N 3 + H.c), (3) 

where H\H 2 = Hi T ia 2 H 2 , a 2 is the second Pauli matrix and Hi T = (Hi° Hi~), 
H 2 T = (H 2 + H 2 °). The parameters m Hi , m N , A\ :k are associated with soft SUSY 
breaking terms.) 2 ] We shall take the parameters A, k, A\ and A k to be real and positive, 
which is a sufficient condition for the vacuum to conserve CP and leads to a choice 



of vacuum in which all the VEVs x, Vi,v 2 are real and positive [Tl]. Note that the 
range of the parameters is restricted by the condition that the vacuum does not break 
QED, which is not automatic in the NMSSM, and is equivalent to the condition that 
m 2 > 0, where m c is the mass of the physical charged Higgs H . Also it is not 
automatic that the vacuum does not break QCD although it has been checked that 
this is the case for values of r and tan (5 in the range 0.05-20 ||11|| . Finally, the range of 



2 In our analysis we shall not consider any other soft SUSY breaking parameters. 



the parameter Ak is restricted by the requirements that < Vm ggs >< 0, and all mass 
squared eigenvalues are positive. We shall use the above conditions in our analysis. 

The effective theory below M SUSY is just the standard model with two light Higgs 
doublets and a light Higgs singlet. Thus the Higgs potential at some low energy scale 
/i < Msusy is given by the general expression 3 

V mgg s = ^A 1 (F 1 t F 1 ) 2 + ^A 2 (^ 2 t F 2 ) 2 + (A 3 + A 4 )(i/i t i/i)(i/2 t ^2) 

- A 4 |# 2 t # 1 | 2 + AslAfl^l 2 + A 6 |Af|# 2 | 2 
+ A 7 (A^ 2 #i# 2 + #.c.) + A 8 |iV| 4 

+ mi 2 \H x \ 2 + m 2 2 \H 2 \ 2 + mz 2 \N\ 2 

- m i {HiH 2 N + H.c.) — -m 5 (N 3 + H.c). (4) 

3 

Comparing Eqs.(3) and (4), the running quartic couplings A, and the mass parameters 
rrii must satisfy the following boundary conditions at M SUS y 

Ai = A 2 = i(# 2 2 + gi 2 ), A 3 = ^{g 2 - gi 2 ) 

A4 = A 2 — 2#2 2 ' A5 = A6 = A 2 , A7 = — Xk, Xs = k 2 , 

mi = m Hl , m 2 = m H2 , m 3 = m N , 

m 4 = XA X , m 5 = kA k . (5) 

At energy scales /1 below Msusy, the values of the quartic couplings may be ob- 
tained by solving the following RG equations which we have derived for the potential 



in Eq.(4) with the aid of the general results in ref. [18 



16tt 2 ^ = 12Ai 2 + 4A 3 2 + 4A 3 A 4 + 2A 4 2 + 2A 5 2 
at 

- Xi(3gi 2 + 9g 2 2 ) + + ^ 2 4 + \g x 2 g 2 2 , 

dX 

16tt 2 -^ = 12A 2 2 + 4A 3 2 + 4A 3 A 4 + 2A 4 2 + 2A 6 2 
at 



3 There is a global U(l) symmetry of the quartic potential which forbids the radiative generation 
of other terms not included in Eq.(4). 



A 2 (3( ?1 2 + 9g 2 2 ) + ^x 4 + ^ 2 4 + 
+ 12h t 2 X 2 - I2h t \ 



<9A 

16ir2 ~df = 2(A 1 + A 2 )(3A 3 + A 4 ) + 4A 3 2 + 2A 4 2 2A 5 A 6 

- A 3 (3^ 2 + 9g 2 2 ) + ^x 4 + ^ 2 4 - ^iV + 6^ 2 A 3 , 
16tt 2 -^ = 2A 4 (Ai + A 2 + 4A 3 + 2A 4 ) + 4A 7 2 

- A 4 (3^ 2 + 9g 2 2 ) + 30i V + 6/i t 2 A 4 , 

<9A 

16tt 2 -^ = 2A 5 (3A! + 2A 5 + 4A 8 ) + 2A 6 (2A 3 + A 4 ) + 8A 7 2 



dt 



2 d\ 6 



^\ 5 (3 gi 2 + 9g 2 2 ), 



~dt = 2As ( 2A3 + A4 ) + 2A 6 (3A 2 + 2A 6 + 4A 8 ) + 8X/ 
- ^A 6 (3( ?1 2 + 9( ?2 2 ) + 6h t 2 X 6 , 



16tt 2 -/ = 2A 7 (A 3 + 2A 4 + 2A 5 + 2A 6 + 2A 8 ) 
at 

- h 7 (3 gi 2 + 9g 2 2 ) + 3h 2 X 7 , 

IGtt 2 ^ = 2A 5 2 + 2A 6 2 + 4A 7 2 + 20A 8 2 , (6) 
at 

where t — In /i, and h t obeys the RG equation 

16vr 2 ^ = \h? - h t (8g, 2 + 9 -g 2 2 + ^ 9l 2 ). (7) 

The gauge couplings gi obey the RG equations of the standard model with two Higgs 
doublets and three fermion families below M S usy, namely 

16tt 2 -^ = -Gig* (ci = -7, c 2 = 3, c 3 = 7), (8) 

and we shall take ai i2 ^(M z ) = 0.0102,0.0336,0.113, respectively. 

The minimisation conditions implied by dv ^ t99a = and 9Vhi " s = allow us to 

L J ovi ox 

eliminate the low energy parameters m 4 , m 2 , m 3 . The remaining masses m 4 and m 5 
are related to the parameters A x and A k at M S ijsy by Eq.(5). Below this scale we 
shall regard m 4 and as free parameters. The charged Higgs squared mass is given 



by 

2x 

m c 2 = ■ o/j ( m 4 - A 7 x) - A 4 w 2 , (9) 
sm2p 

which allows us to eliminate the parameter m 4 in favour of the more physical param- 
eter m 2 . The physical neutral CP-odd pseudoscalar mass squared symmetric matrix 
elements are 

v 2 3 
(M p 2 J n = 3xm 5 + ^ sin 2 2/3(m 2 + X 4 v 2 ) - -X 7 v 2 sin 2(3, 

(M 2 J 12 = |-sin2/3(m 2 + A 4 t; 2 ) + 3t;a;A 7 , 

(M 2 J 22 = m 2 + A 4 t; 2 . (10) 

The neutral CP-even scalar mass squared symmetric matrix elements are, in the basis 
1,2,3 = H^H^N, 

(M 2 ) n = 2AiV + (m 2 c + X 4 v 2 )sm 2 (3, 

(M 2 ) 12 = sm2(3[(\ 3 + ±\ 4 )v 2 -±m 2 c ], 

(M 2 ) 22 = 2A 2 v 2 2 + cos 2 / 3(m 2 + A 4 t; 2 ), 

(M 2 ) 13 = x(2A 5 wi + A 7 t; 2 )-^|sin2/3(m 2 + A 4 t; 2 ), 

(M 2 ) 23 = x(2A 6 v 2 + A 7 t;i)-|^sin2/3(m 2 + A 4 t; 2 ), 

v 2 1 
(M 2 ) 33 = 4A 8 x 2 + — sin 2 2/3(m 2 + A 4 t; 2 ) + -A 7 v 2 sin2/5-xm 5 . (11) 

The model is then specified by the 6 parameters A, k, tan/5 = v 2 /vi, r = x/v, m 5 , 
and m 2 . Our detailed procedure for evaluating Higgs masses is as follows. Given 
the gauge couplings gi(Mz), we find gi{m t ) by solving the RG equations (8). For a 
given m t and tan/3 we obtain h t {m t ), then use its RG equation (7) and those for the 
gauge couplings (8) to evaluate h t (M SUSY ) and gi(M SUSY )- Given these couplings 
and some choice of A and k we then use the boundary conditions in Eq.(5) and the 
RG equations (6), (7), (8), to find Aj(/i), where /i is close to mass-shell for the Higgs 
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boson of interest. [] In practice we shall take Msusy = 1 TeV, /i = 150 GeV. The 
neutral Higgs masses are then calculated from Eqs.(10),(ll) with Aj = Ai(yu). The 
diagonalisation of these matrices also gives the mixing angles, so that the amplitude 
of the singlet state N in each of the neutral mass eigenstates may be calculated. 

We shall now obtain an upper bound on the mass of the lightest neutral scalar 
Higgs boson h°, by using the fact that m\ must not exceed the lower eigenvalue of 
the upper 2x2 block matrix, (M 2 )^-, i,j = 1,2. The resulting upper bound is a 
complicated function of m 2 

m h 2 <^(A + m 2 ) l -^{ m( ? + BY + C^-B\ (12) 

where 

A = v 2 (\ 1 + \ 2 + \ 4 ) + v 2 (\ 1 - \ 2 )cos2(3, 

B = -w 2 cos2/5[(Ai-A 2 ) + (Ai + A2-A 4 )cos2/3]-2t; 2 sin 2 2/5(A3 + ^A4), 
C 2 = v 4 [(\ 1 -\ 2 ) + (\ 1 + \ 2 -\ A )cos2(3] 2 + 4v 4 sm 2 2(3(\ 3 + ^\ i ) 2 . (13) 

It is easy to show that the above bound reaches a maximum asymptotically for m c — > 
oo. Since {C 2 — B 2 ) > 0, we obtain the m c -independent bound 

m h 2 <^(A-B). (14) 

Inserting A, B from Eq.(13) and using the boundary conditions in Eq.(5) leads to an 
upper bound of the form 

m h 2 < M z 2 + (AV - M z 2 ) sin 2 2(3 
v 2 

+ y + sx 2 + 25A 3 + 25A 4 ) + 2(5\ 1 - 5X 2 ) cos 2(3 

+ (5A 1 + 5A 2 -25A 3 -25A 4 )cos 2 2/3], (15) 

4 In integrating the RG equations from Msusy down to /x we have decoupled the terms involving 
ht below the scale m t , but we have not decoupled any Higgs bosons below their mass scales. We 
would not expect such a decoupling to change our results significantly providing there are no very 
heavy scalars. 
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where 5Xi = Aj(/i) — Xi(MsusY), and, strictly speaking, Mz should be evaluated at 
Msusy- Using the RG equations (6), keeping only terms quadrilinear in A, h t , and k 
and making a small <5Aj approximation, Eq.(15) yields the analytic bound f] 



m h 2 < M z 2 + (AV - M z 2 ) sin 2 2(3 

327T 2 I /J 



+ ^ ln ( ) K 12 ^ 4 - 12A V - 8A2fc2 - 24A ") - 24/i * 4 cos ^ 



+ (12V + 12A 2 V + 8AV + 8A 4 )cos'2/3]. (16) 

In Eq.(16), the tree level bound in Eq.(2) is corrected by a function of h t , cos 2/3, A 
and k. It is easy to show that this function is maximised for A > X max and k = 0, 
and so we shall use X ma x and k = 0, as in the case of the tree level bound discussed 
previously. With this information in hand, let us now return to the more exact result 
in Eq.(15) in which the values of 5Xi are obtained from numerically integrating the RG 
equations, decoupling the top quark below its mass and choosing fi = 150 GeV. For 
each value of m t , we have determined the value of h t (M SUSY ) and the corresponding 
value of X max which maximise the r.h.s. of Eq.(15) from a numerical analysis of the 
triviality condition as discussed previously. These are shown in Fig.l, together with 
sin (3 which is eliminated from Eq.(15) using m t = h t (m t )v sin/?. 

Using the parameters in Fig.l, we obtain the upper bound on the lightest neutral 
Higgs scalar in the NMSSM as a function of m t as shown in Fig.2. We again emphasise 
that we have considered the more general case of two light Higgs doublets plus a 



light singlet all below the scale Msusy-, whereas previous results [16, 17| have only 



considered the case of one light Higgs boson. From Fig.2 we find that for m t = 90 



GeV, rrih < 145 GeV, which agrees to within 2 GeV with both refs.|L6], [I7|. Our 
bound for larger top quark masses, although qualitatively similar to refs. [jl6| , p~7f| , can 
differ by several GeV. It can easily be checked that the analytic formula in Eq.(16) 
agrees with the full numerical bound shown to within 5 GeV across the entire m t 



5 Note that Eq.(16) reduces to the bound in the MSSM in the limit that A -> 0. 
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region, giving a bound somewhat lower than the true bound for m t = 90 — 160 GeV, 
somewhat higher for m t = 180 — 190 GeV, and in good agreement for m t ~ 170 GeV. 
If the A dependent terms in the radiative corrections in Eq.(16) were dropped, this 
would produce a bound about 7 GeV higher than the true bound over the entire m t 
range. 

We shall now discuss the physical parameter regions in which the upper bound 
on rrih in Fig.2 is approximately realised, and discuss the light Higgs boson spectrum 
in these regions. As mentioned earlier, the m c independent bound in Eqs.(15),(16) 
is exactly valid only asymptotically, but actually the bound rapidly approaches its 
maximum for large values of m c of the order of the bound itself. Since it is known that 
the allowed range of m c is controlled by r and that larger values of m c are associated 
with the larger values of r, we shall restrict ourselves to r = 1.0 and 10.0, as shown 
in Figs. 3 and 4. 

Fig. 3 (a) shows the light Higgs boson spectrum as a function of m c for m t = 150 
GeV, for r = 1.0, where tan/3 = 1.7 is from Fig. 1. We have taken A = 0.65 and 
k = 0.1, rather than \ ma x = 0.7 and k = since these values would imply a massless 
CP-odd scalar, which is phenomenologically ruled out. We have also taken m 5 = 0. 
Only the lightest neutral CP-even scalars (solid) and the lightest CP-odd pseudoscalar 
(dashed) are included in Fig. 3 (a). The lightest scalar mass rrih reaches a maximum of 
m h = 88 GeV. In Fig. 3(b) we show the N components of the bosons in Fig. 3 (a). This 
is of interest since the N component is a gauge singlet, and thus decoupled from the 
gauge bosons. The lightest CP-even scalar is about 80% decoupled rising to being 
almost 100% decoupled at its maximum. It is obvious that the bound which follows 
from the upper 2x2 block matrix can only approach saturation if the full scalar matrix 
is approximately block diagonal. Since the lightest scalar is approximately decoupled 
near its maximum mass, the physically observable second lightest scalar must respect 
the bound at this point. Taking m 5 to be non-zero and positive (negative) serves to 
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lower (raise) (M s 2 ) 33 (identified as the decoupled scalar mass) while leaving the other 
elements unchanged, as is clear from Eq.(ll). Finally it is clear from Figs. 3(a), 3(b) 
that the lightest CP-odd state presents no phenomenological problems, since as well 
as being heavy it is over 90% decoupled over the entire m c range. 

In Fig. 4(a) we show the light Higgs spectrum for r = 10.0, choosing the other 
parameters to be the same as in Fig. 3(a). The relevant range of m c is now 2.3-2.5 
TeV which is strictly outside the range of validity of our calculation, since we do 
not consider Higgs decoupling and choose \i = 150 GeV (of order rrih) as in Fig.3(a). 
Nevertheless Fig.4(a) shows that the lightest CP-even scalar mass reaches a maximum 
of 118 GeV, close to the minimum of the physically observable second lightest scalar 
in Fig. 3 (a). Fig.4(b) shows that in this case the lightest scalar is the detectable 
one, being almost 0% decoupled at its maximum, while the second lightest scalar is 
almost 100% decoupled at this point. The CP-odd boson in Fig. 4 is now nearly 100% 
decoupled over the whole range of m c . 

In conclusion, the lightest CP-even neutral Higgs boson in the NMSSM must 
respect the bounds, rrih < 145 GeV for m t = 90 GeV, decreasing to rrih < 123 
GeV for m t = 180 GeV. Our analysis is based on the assumption of a SUSY desert 
between M S usy = 1 TeV and Mqut = 10 16 GeV, and involves the assumption of hard 
decoupling of the superpartners. However, our calculations do include the effects of 
other light Higgs bosons which previous analyses have ignored. We have presented 
an analytic expression in Eq.(16) for the radiative corrections in the NMSSM, which 
agrees well with our full numerical bound in Fig.2. In addition we have discussed 
the light Higgs boson spectrum in a region of parameter space which approximately 
realises the bound, and have seen that there are indeed other light Higgs bosons whose 
effects must be taken into account. For r ~ 1 (Fig.3), the lightest CP-even scalar is 
not very strongly physically coupled when it is close to its maximum mass, and the 
second lightest CP-even scalar, which must also respect the bound, may be easier to 
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discover. 
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Figure Captions 



Figure 1 : The values of h t (solid), \ max (short dashed) and sin/? (long dashed) 
realising the upper bound on rrih in Figure 2. 

Figure 2 : The full numerical bound on the lightest neutral CP-even Higgs scalar 
as a function of m t . 

Figure 3(a) : The light Higgs boson spectrum for the parameters m t =150 GeV, 
r = 1.0, tan/? = 1.7, A = 0.65, k = 0.1. The solid lines are CP-even scalar masses, 
and the dashed line is the CP-odd mass. 

Figure 3(b) : The iV components of each of the bosons in Fig.3(a) where the 
solid line is for the lightest CP-even scalar, the dashed line is for the second lightest 
CP-even scalar, and the dotted line is for the CP-odd pseudoscalar. 

Figure 4(a) : The light Higgs boson spectrum for r = 10.0. The other parameters 
are as in Figure 3(a). The solid lines are CP-even scalar masses, and the dashed line 
is the CP-odd mass. 

Figure 4(b) : The iV components of each of the bosons in Fig. 4(a) where, as in 
Fig.3(b), the solid line is for the lightest CP-even scalar, the dashed line is for the 
second lightest CP-even scalar, and the dotted line is for the CP-odd pseudoscalar. 
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